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Study guide 


There are five sections in this chapter. They are 
intended to be studied consecutively in five study 
sessions. 


Section 1, which starts and ends with video 
sequences, requires the use of a DVD player, and 
Section 4 requires the use of the computer together 
with Computer Book B. 


The pattern of study for each session might be as 
follows. 


Study session 1: Section 1. 
Study session 2: Section 2. 
Study session 3: Section 3. 
Study session 4: Section 4. 
Study session 5: Section 5. 


Each of the sections should take two to three hours 
of study. 


Section 1 introduces matrices, and Section 2 covers 
a little of their theory. Sections 3 and 4 contain a 
modelling case-study that shows one application of 
matrices, whilst Section 5 continues with more 

of the theory of matrices. An alternative study 
pattern, which you might prefer, is to study 
Section 5 before Section 3, thus putting all the 
theory together before studying an application. 


In addition to the work outlined above, there is an 
optional band on CDA5505, in which various 
applications of matrices are described. The best 
time to listen to this band is after you have 
completed your study of the chapter. 


The optional Video Band B(iii) 

Algebra workout — Arithmetic with matrices 
could be viewed at any stage during your 
study of this chapter. 


Introduction 


As its title suggests, this chapter introduces and uses mathematical objects 
called matrices. Matrices are a very important topic in mathematics and 
find widespread use in mathematical modelling. In fact, almost every 
problem that is too complicated to be solved with paper and pencil is first 
reduced to a matrix problem before being solved numerically by computer. 


Matrices are introduced in Section 1 through the concept of networks. The 
network image provides a way of visualising matrices and how to multiply 
them. In Section 2, the matrix concept is built upon and standard matrix 
notation is defined. 


In Section 3, a population problem is studied. This problem can be solved 
using matrices with the basic operations defined in Section 2. The model 

which we shall use is a sequence model similar to one of those used in the 
previous chapter, but with the population split into subpopulations. 


Section 4 is a computer session designed to show how useful computer 
software is for dealing with matrices. The section starts by introducing the 
ways in which the computer can manipulate matrices, then proceeds to 
investigate further the population problem from Section 3. 


In Section 5, pairs of simultaneous linear equations are expressed in terms 
of matrices and solved, where possible, by manipulating those matrices. 


The singular of matrices is 
matria. 


1 Networks and matrices 


Where no confusion arises, 
the word ‘network’ is used 
to mean either ‘network 
diagram’ or ‘physical 
network’. Similarly, the word 
‘pipe’ will mean either ‘pipe 
in a network diagram’ or 
‘physical pipe’. 


Note that two pipes that 
cross on a diagram do not 
join with each other unless 
the crossing point is marked 
as a node, by means of a 
circle. For example, the pipe 
from node A to node V does 
not join with the pipe from 
node B to node U, because 
the crossing point is not 
marked as a node. 


To study this section you will need a DVD player and DVD00104. The 
section starts with a video sequence (Subsection 1.1) and concludes with 
one (Subsection 1.3). 


In this section you will meet mathematical objects called matrices, and see 
how to ‘multiply’ two matrices. Matrices can be used to represent 
particular kinds of network, and ‘combining’ two networks corresponds to 
matrix multiplication. 


1.1 Networks 


This section starts with a discussion of networks, on the video band. 


Watch video band B(ii), ‘Networks and matrices’, until you reach the 
caption ‘Stop the DVD’. 


Two of the main ideas from the video are the following. 


© A network diagram is a mathematical representation of a physical 
network. Each point at which a network diagram branches is called a 
node. Two nodes in a network can be connected by a pipe. 


© In the networks that we consider, the pipes are labelled by numbers. 
Suppose that 1 litre of water is input at a node, and that it flows from 
the node down a number of pipes. The number labelling each pipe is 
the amount (in litres) of that 1 litre of water which emerges via that 
pipe at the node at the bottom of the pipe. 


The first network introduced in the video can be represented by the 
following network diagram. 


Figure 1.1 Network with two input nodes (A and B), three output nodes 
(U, V and W) and six pipes 


Consider what happens when 1 litre of water is input at the node 

labelled A. The water flows downwards, and can be collected at the nodes 
labelled U, V and W. The amount of water reaching each node can be 
read off the network diagram by using the pipe labels. So in this case 

0.4 litres of water reach node U, 0.2 litres of water reach node V, 

and 0.4 litres of water reach node W. 


Now try the following activity. 


SECTION 1 NETWORKS AND MATRICES 


Activity 1.1 Calculating network outputs 


Consider the network diagram shown in Figure 1.1. 


(a) How much water is output at the nodes U, V and W if the only input 
is 1 litre of water at node B? 


(b) How much water is output at the nodes U, V and W if the only input 
is 2 litres of water at node A? 


Solutions are given on page 46. 


In Activity 1.1(b), you should have found that doubling the amount of Activity 1.1(b) indicates that 
water input doubles the amount of water at each output node. The same the label on the pipe from 
principle applies if water is input at both node A and node B: the amount node A to node U (0.4), for 
of water in an output node is the sum of the water reaching that node example, gives the proportion 
from node A and from node B. Consider, for example, what happens when © water input at node A that 
1 litre of water is input at node A and another 1 litre is input at node B. reaches node U. The other 
At node U, 0.4 litres of water arrive from node A and 0.3 litres of water pipe Telbele) Gam De Iverereted 
arrive from node B. So, in total, there will be 0.7 litres of water in the Me Nera Oo 

container at node U. Similarly, 0.2 + 0.1 = 0.3 litres of water collect at 

node V, and 0.4 + 0.6 = 1 litre of water collects at node W. 


Now try using this reasoning yourself to calculate network outputs. 


Activity 1.2 Calculating network outputs 


Consider the network diagram shown in Figure 1.1. 


(a) How much water is output at the nodes U, V and W if 1 litre of water 
is input at node A and 2 litres of water are input at node B? 


(b) How much water is output at the nodes U, V and W if x litres of 
water are input at node A and y litres of water are input at node B? 


Solutions are given on page 46. 


Activity 1.3 Calculating more network outputs 


Consider the following network diagram, which represents the second 
network shown in the video. 


Figure 1.2 Network diagram with three input nodes and two output nodes 


(a) How much water is output at nodes C and D if 1 litre of water is 
input at node U and 2 litres of water are input at node V? (There is 
no input at node W.) 


(b) How much water is output at nodes C and D if x litres of water are 
input at node U, y litres of water are input at node V, and z litres of 
water are input at node W ? 


Solutions are given on page 46. 


CHAPTER B2 MODELLING WITH MATRICES 


At the end of the video sequence you saw the effect of putting together the 
two networks represented in Figures 1.1 and 1.2. This combined network is 
shown in Figure 1.3. 


Figure 1.3 The combined network shown in the video 


Using this diagram of the combined network, the amount of water flowing 
from the input nodes to the output nodes can be calculated. For example, 
suppose that 1 litre of water is input at node A (and there are no other 
inputs). Then we know from the discussion preceding Activity 1.1 that the 
amounts of water passing through nodes U, V and W are 0.4, 0.2 and 0.4 
litres, respectively. These ‘outputs’ from the top network can be viewed as 
‘inputs’ to the bottom network. To show how this is done, consider the 
output node C’.. From the network diagram it can be seen that: 


© 0.2 of the water passing through node U reaches node C, 
so 0.2 x 0.4 = 0.08 litres of water from node A reach node C 
via node U; 


© 0.6 of the water passing through node V reaches node C, 
so 0.6 x 0.2 = 0.12 litres of water from node A reach node C 
via node V; 


© 0.3 of the water passing through node W reaches node C, 
so 0.3 x 0.4 = 0.12 litres of water from node A reach node C 
via node W. 


To find the total amount of water output at node C, we add the 
contributions from the three different routes from A, to obtain 
0.08 + 0.12 + 0.12 = 0.32 litres of water output at node C. 


Activity 1.4 Flow through a combined network 


Both parts of this activity refer to the combined network as shown in 
Figure 1.3. 


(a) Calculate the amount of water output at node D if the only input is 
1 litre of water at node A. 


(b) Calculate the amount of water output at nodes C and D if the only 
input is 1 litre of water at node B. 


Solutions are given on page 46. 


SECTION 1 NETWORKS AND MATRICES 


The combined network can be considered as equivalent to a simpler 
network with two inputs at nodes A and B and two outputs at nodes 
Cand D. The above calculations would be a lot simpler if we could find 
the network diagram for this simpler network. In fact, that is what we 
have just done. Remember that in the video, when numbers were first 
associated with pipes, the number that labelled a pipe was the amount of 
water that flowed down the pipe when 1 litre of water was input to the 
node at the top. 


From the calculations preceding Activity 1.4, if 1 litre of water is input at 
node A, then 0.32 litres of water will be output at node C. This means that 
the pipe from node A to node C in the simpler network can be labelled by 
the number 0.32. The calculations in Activity 1.4 provide the labels for the 
other three pipes, giving the network diagram shown in Figure 1.4. 


Cc D 


Figure 1.4 A two-input, two-output network which is equivalent to the combined 
network shown in Figure 1.3 


To see how much easier the calculations are with this new network, try the 
following activity. 


Activity 1.5 Flow through a two-input, two-output network 


All parts of this activity refer to the two-input, two-output network shown 
in Figure 1.4. 


(a) Calculate the amount of water output at nodes C' and D if the only 
input is 1 litre of water at node B. 


(b) Calculate the amount of water output at nodes C' and D if « litres of 
water are input at node A and y litres of water are input at node B. 


Solutions are given on page 47. 


1.2 From networks to matrices 


In Subsection 1.1 you calculated the amount of water leaving a pipe 
network when various quantities of water were input. The calculations 
were lengthy to write out because you constantly needed to write phrases 
such as ‘input at node A’ or ‘output at node U’ or ‘amount reaching 
node C' via node U’. Mathematicians generally try to avoid superfluous 
words, so these calculations are usually expressed more succinctly using 
matrices, which are introduced in this subsection. 


Figure 1.5 The first network 
seen in the video 


Vectors will be used to 
represent other quantities of 
interest in later sections. 


Consequently, we shall refer 
to ‘input’ and ‘output’ 
vectors. 
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We start by considering again our first network, shown in Figure 1.5. This 
network has two inputs and three outputs. The amount of water input at 
the two input nodes can be represented by a column of two numbers. For 


example, the column & represents an input of 1 litre of water at node A 


0 
and 0 litres of water at node B. So we have replaced the phrase ‘litres of 
water input at node A’ by a specified position in a column of numbers by 
ordering the nodes (in this case alphabetical order). The column of 
numbers above is called a vector, and the numbers appearing in the 
column are called its components. In the above example, | is called the 
first component and 0 is called the second component. 


Here are some more examples of vectors. 


0 1 x o 
GROROM(: 
0.4 
In the context of the network in Figure 1.5, the first of these vectors 
represents 0 litres of water input at node A and 1 litre of water input at 
node B. The second vector represents 1 litre of water input at node A and 
2 litres of water input at node B. The third vector is the general case 
of x litres of water input at node A and y litres of water input at node B. 
The fourth vector could be used to represent the outputs from the network 
—namely, 0.4 litres of water output at node U, 0.2 litres of water output at 


node V, and 0.4 litres of water output at node W — corresponding to some 
vector of inputs. 


Try the following activity, which reinforces the link between vectors and 
the inputs to and outputs from networks. 


Activity 1.6 Representing inputs and outputs by vectors 


All parts of this activity refer to the first network that we considered, as 
shown in Figure 1.5. 


(a) How can inputs of 2 litres of water at node A and 3 litres of water at 
node B be represented by a vector? 


(b) What is the first component of the vector that you wrote down as your 
answer to part (a)? 
0.3 
(c) What outputs from the network are represented by the vector | 0.1 |? 
0.6 
Solutions are given on page 47. 


Having seen that the inputs to and outputs from a network can be 
represented by vectors, we turn our attention to representing the network 
itself. How do we represent the way in which the outputs are related to the 
inputs? The answer is to use a matrix, so we now proceed to define this 
concept. 


SECTION 1 NETWORKS AND MATRICES 


A matrix is a rectangular array of numbers, which can be used to hold The word ‘matrix’ is usually 
quantitative information in a structured way. Examples are pronounced ‘may-trix’. 

1 2 -1 aA 1 2 

5 3 0 3.4)" 
Each number in a matrix is called an element of that matrix. So, for Note that although one 
example, the first matrix above has six elements and the second matrix has sometimes talks about the 
four elements. elements of a vector, one 


never talks of the components 


The elements of a matrix are arranged in rows and columns. A row of a oh Ate 


matrix is a horizontal line of numbers in the matrix. For example, 
(1 2 -—1) is the first row of the first matrix written above, and 
(5 3 0) is the second row of that matrix. A column of a matrix is a 
vertical line of numbers in the matrix. For example, the first matrix 

: ‘ 1 2 —1 
written above has three columns, which are ( 5 ), ( | and ( 0 ). 
There is a similarity between the columns of a matrix as shown above and 
the vectors defined earlier: this is no coincidence. In this chapter, ‘vector’ 
is simply a shorter name for ‘matrix with one column’. In the next 
chapter, you will meet vectors in a different context. 


In order to represent a network by a matrix, the numbers associated with 
each pipe should appear in the matrix. There are several ways of doing 
this. By convention, the network in Figure 1.5 is represented by the matrix 


0.4 0.3 
0.2 0.1 
0.4 0.6 


You can see that the elements in a column specify the flow from a node 
(A or B), and the elements in a row specify the flow to a node (U, V 
or W). The reason for adopting this convention will be apparent in the 
discussion following Activity 1.7. 


Now see if you can apply this convention to represent the second network 
that you saw in the video. 


Activity 1.7 Representing a network 


The second network that you saw in the video is shown in Figure 1.6. 


(a) Write down a vector which represents an input of 1 litre of water at 
node V, and no other inputs. 


(b) How many components does the input vector have? How many 
columns do you expect the matrix representing the network to have? 


(c) How many components does the output vector have? How many rows 


do you expect the matrix representing the network to have? Figure 1.6 The second 


(d) Using your answers to parts (b) and (c), write down the matrix which network seen in the video 
represents this network. 


Solutions are given on page 47. 


Th 


The following discussion 
relates to the network in 
Figure 1.5, which was 
represented by the matrix 


0.4 0.3 
0.2 0.1 
0.4 0.6 


IZ 
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Now we look at some of the results calculated earlier for networks, and see 
how they translate into matrix notation. The first thing we did was to 
calculate the outputs from the network in Figure 1.5 when 1 litre of water 
was input at node A. The inputs are represented by the vector 


(0). 


The outputs calculated were 0.4 litres of water at node U, 0.2 litres of 
water at node V, and 0.4 litres of water at node W, which can be 
expressed in vector form as 


0.4 
0.2 
0.4 


How do we represent the relationship between the input vector and the 
output vector, using matrices? We take our cue from the way that 
functions are represented. For example, if a function f takes an input 
number, x, to an output number, y, then this can be written as f(x) = y. 
By analogy, we represent the relationship between the input vector and the 
output vector as 


0.4 0.3 1 0.4 

0.2 0.1 ( 0 ) = | 0.2 

0.4 0.6 0.4 
The interpretation of this equation is that the input vector in the middle is 
acted upon by the matrix to give the output vector on the right. 


The above equation was derived from results calculated using network 
diagrams. What is the rule for calculating the output vector given the 
matrix of a network and an input vector? This is what we now proceed to 
find. 


In Activity 1.2 you considered the general case of x litres of water input at 
node A and y litres of water input at node B, and found that 


0.42 + 0.3y litres of water are output at node U, 

0.2% + 0.1y litres of water are output at node V, 

0.4% + 0.6y litres of water are output at node W. 
This corresponds to the output vector 


0.42 + 0.3y 
0.27 + 0.ly 
0.4a + 0.6y 


So, in matrix notation, the relationship between the input vector and the 
output vector is 


0.4 0.3\ / 0.42 + 0.3y 
0.2 0.1 ( ) = | 0.2 +0.ly 
0.4 06/ \¥ 0.42 + 0.6y 


Now look at this equation. The first component of the output vector is the 
sum of the products of each of the numbers in the first row of the matrix 
with the corresponding component in the input vector, as highlighted in 
Figure 1.7. 


SECTION 1 NETWORKS AND MATRICES 


0.4 0.3 ” O.4e + O.3y 
0.2 0.1 ( ) = | 0.244 0.ly 
0.4 0.6 0.42 + 0.6y 


Figure 1.7 Obtaining the first component of the output vector 


Similarly, the second and third components of the output vector are 
obtained from the second and third rows of the matrix, respectively, as 
shown in Figure 1.8. 


0.4 0.3 ” O.4e + O.3y 0.4 0.3 ‘i O.4e 4+ O.3y 
0.2 0.1 ( ) = | 0.244 0.1ly 0.2 0.1 ( ) = | 0.244 0.ly 
0.4 06/ \¥ 0.42 + 0.6y 0.4 0.6 0.42 + 0.6y 


Figure 1.8 Obtaining the second and third components of the output vector 


So we can find any component of the output vector from the elements in 
the corresponding row of the matrix and the components of the input 
vector, by forming the sum of their products. This operation is called 
matrix—vector multiplication. Note that a multiplication sign is not 
used to denote this operation; the matrix and vector are simply written 
next to each other. 


This subsection ends with activities which ask you to find the output from 
a network when given the matrix representing the network, instead of the 
network diagram. The first activity refers to a familiar network for which 
you have the network diagram to guide you; the second activity refers to a 
new network. 


Activity 1.8 Using matrices to calculate outputs 


The second network that you saw in the video is shown in Figure 1.6. 


(a) Write down a vector which represents an input of 1 litre of water at 
node U, and no other inputs. 


(b) The matrix representing the network was found in Activity 1.7. Use 
this matrix to calculate the output vector for the above input vector, 
and interpret your answer in terms of the amounts of water collected 
at nodes C' and D. 


Solutions are given on page 47. 


Activity 1.9 Finding outputs 
A network has two input nodes and two output nodes. Its matrix is 
0.33 0.25 
0.67 0.75) ° 
Find the output vector for each of the following input vectors. 


©() o(8) © (8) 


Solutions are given on page 47. 
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This expression for the output 
must evaluate as a vector 
with two components, but we 
do not yet know how the two 
matrices on the left combine. 
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Matrices were invented by Arthur Cayley (1821-95). In 1858, in a 
paper presented to the Royal Society, he gave the rules for matrix 
operations (covered in Section 2) and the conditions under which a 
matrix has an inverse (covered in Section 5). However, it was Cayley’s 
friend and fellow mathematician James Joseph Sylvester (1814-97) 
who coined the term ‘matrix’, in 1850. 


1.3 From combining networks to matrix 
multiplication 


In Subsection 1.1 you saw that two networks can be combined to form a 
simpler equivalent network. To be more precise, you saw how the first two 
networks on the video can be combined to produce the two-input, 
two-output network shown in Figure 1.9. 


Figure 1.9 Two equivalent networks 


It turns out that this simplification can also be achieved by ‘combining’ 
matrices. Suppose that the input to the left-hand network shown in 
Figure 1.9 is the vector i Then the output from the top part of the 


network can be written in the form 


U 0.4 0.3 - 
oe )=102 01 ( ). (1.1) 
w 0.4 06/ \¥ 


This three-component output vector is used as the input to the bottom 
part of the left-hand network. Thus the output from the network is given 
by the following matrix—vector multiplication: 


0.2 0.6 0.3) ( ° 
0.8 0.4 0.7 , 
w 
We can substitute for the three-component vector using equation (1.1), to 
obtain the output from the whole network in terms of the input: 


(3 0.6 a) eas (2) Ae) 
0.8 0.4 0.7 04 06 y 

Now consider the right-hand network shown in Figure 1.9, which is 
represented by the matrix 


0.382 0.3 
0.68 0.7 /)° 


SECTION 1 NETWORKS AND MATRICES 


The output from the network for the input (5) is 


(oss or) (7): (13) 


Since the two networks are equivalent, expression (1.2) must be equal to 
expression (1.3). This will happen if there is a sense in which meaning can 
be given to the left-hand side of the matrix equation 


0.2 0.6 0.3 e ne _ {0.32 0.3 LA The notion of equality 
0.8 0.4 0.7 . i ~ \ 0.68 0.7 /° (1.4) between matrices is defined in 
0.4 0.6 Section 2. 


Equation (1.4) was deduced by considering the combination and 
simplification of two networks. We now define a method of ‘combining’ 
matrices that is consistent with equation (1.4). This method of 
‘combining’ matrices is called matrix multiplication, and the ‘combined’ 
matrix is called the product matriz. 


We start by considering the element 0.32 in the top left-hand corner of the 


soe a How was this calculated? If you look back 
to Subsection 1.1 (page 8), then you will see that we added the 
contributions flowing down each pair of pipes that could lead from node A 
to node C.. This is highlighted in the left-hand diagram of Figure 1.10. In 
the right-hand diagram of Figure 1.10, the corresponding matrix elements 
are highlighted. 


product matrix ( 


(= 0.6 a) Fe A 
0.8 O04 O.7 0.4 06 


Figure 1.10 Highlighting shows pipes and corresponding matrix elements 


Looking at the right-hand diagram in Figure 1.10, we see that the top left 
entry of the product matrix is formed from the top row of the first matrix 
and the left-hand column of the second matrix. What is the rule for 
combining a row and a column? The answer is exactly the same as when 
calculating a matrix-vector product: add together the products of 
corresponding elements. In this case, we obtain 


(0.2 x 0.4) + (0.6 x 0.2) + (0.3 x 0.4) = 0.32. 


Now let us see this idea at work again. How is the bottom left element, 
0.68, of the product matrix formed? In terms of combining networks, it is 
the sum of the contributions flowing from node A to node D via any of the 
three possible routes. In terms of the matrix product, we expect that the 
bottom left element of the product matrix should be formed from the 
bottom row of the first matrix and the left-hand column of the second 
matrix. The network routes and matrix elements are highlighted in 

Figure 1.11, overleaf. 
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Matrix manipulation is the 
subject of Section 2. 
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e 0.6 i ae ie 
0.8 O04 O07 0.4 0.6 


Figure 1.11 Highlighting shows pipes and corresponding matrix elements 


From the elements highlighted in the matrix multiplication above, we can 
calculate the corresponding element of the matrix product: 


(0.8 x 0.4) + (0.4 x 0.2) + (0.7 x 0.4) = 0.68. 


The process of combining two matrices as described above is known as 
matrix multiplication. The result of multiplying two matrices is called 
the product matrix. 


The next activity asks you to verify that the other two elements in the 
0.32 0.3 


0.68 are correct. 


product matrix ( 


Activity 1.10 Combining networks and multiplying matrices 


(a) Consider the top right element, 0.3, in the above product matrix. 


(i) Sketch the network diagram, and highlight the routes through the 
network that contribute to this element of the product matrix. 


(ii) Write out the matrix multiplication, and highlight the elements of 
the matrices that contribute to this element. 


(iii) Use the highlighted elements of the matrices to re-calculate the 
element of the product matrix. 


(b) Repeat part (a) for the the bottom right element of the product 
matrix. 


Solutions are given on page 48. 


This section ends with another look at finding outputs using matrices and 
multiplying matrices on the video. 


Now watch the video to the end of band B/(ii). 


SECTION 1 NETWORKS AND MATRICES 


Summary of Section 1 


This section has introduced: 


© the idea of a network, and how to calculate the output from a network 
and from combined networks; 


© vectors, as a notation for inputs to and outputs from a network; 


© 


matrices, as a notation for representing networks; 


© matrix—vector multiplication, as a way of calculating the output from 
a network; 


© matrix multiplication, by considering the combination of networks. 


Exercise for Section 1 


Exercise 1.1 


(a) Find the output vector when each of the following input vectors is 
0.35 0.85 
0.65 0.15 /* 


0 (2) (i) (8) 


(b) Sketch the network corresponding to the matrix in part (a), including 
the appropriate label for each pipe. 


acted on by the matrix ( 
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2 Matrix manipulation 


‘2 x 3’ is read as ‘two by 
three’. 


Some people prefer to say 
shape or order rather than 
size here. 
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This section aims to help strengthen your understanding of matrix 
multiplication and introduces some other ways of manipulating matrices. 


2.1 Matrix multiplication 


Section 1 introduced matrix multiplication and in particular evaluated the 
matrix product 


t 0.6 i) 2 oa Ser 7 (2.1) 
0.8 0.4 0.7 0.4 06 0.68 0.7 


There is something special about the numbers of rows and columns of the 
matrices appearing in equation (2.1). In order to discuss this, we need 
some notation. The left-hand matrix has 2 rows and 3 columns. We often 
abbreviate this and say that the left-hand matrix is a 2 x 3 matrix. The 
next matrix has 3 rows and 2 columns, so it is a 3 x 2 matrix. 


In general, an m X n matrix is one with m rows and n columns. We say 
that the size of the matrix is m x n. Also, two matrices of the same size 
are equal if all their corresponding elements agree. 


The special nature of the sizes of the matrices in the above matrix product 
is best seen by writing the sizes next to each other as follows: 


2x3 3x2 gives 2x 2. 


The first observation is that the two 3s are adjacent to each other. There 
is a good reason for this: to combine a row and a column, they must each 
have the same number of elements. The first matrix has 3 elements in each 
row and the second matrix has 3 elements in each column. 


We can indicate the matching of the 3s by drawing a box around them: 


2x/|3 3/x2 = gives 2x 2. 


If we now ignore the numbers in the box, then on the left we have the 
numbers 2 and 2 which appear in the size of the product matrix, 2 x 2. 


In fact, this always happens. Two matrices can be multiplied only if the 
adjacent numbers match. Furthermore, the size of the product matrix is 
given by the remaining numbers when we ignore the adjacent numbers. 


Thus we have a useful method for deciding whether two matrices can be 
multiplied and for obtaining the size of a product matrix. Try the 
following activity to practise this method. 


Activity 2.1 Determining the size of matrix products 


For each of the following, state whether the two matrices can be 
multiplied. If they can, then determine the size of the product matrix. 
(Do not attempt to calculate the product matrix.) 
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3 9 9 
a5 4 4 
(a) 7 3 (oy cl 2 Bis 
=§ 8 ( =o 3) (8) 
i 9 0 10 9 
(c) /-- 43 (d) (3) 2 3) 
—6 8 =o 4 7 


oc s(t 


Solutions are given on page 48. 


Before continuing, here is the general definition of matrix multiplication, 
which summarises what we have found above. 


Matrix multiplication 


Let A and B be two matrices. We can form the product matrix AB Note that we use the 

only if the number of columns of A equals the number of rows of B. convention of upper-case bold 
letters to denote matrices. 
When writing these letters by 
hand, they are usually 


The product matrix AB is formed by ‘combining’ rows of matrix A 
with columns of matrix B, where ‘combining’ means adding the 


products of corresponding elements. The rule for where to put the written as plain upper-case 
resulting numbers is that the element in the 7th row and jth column letters (with no underlining 
of AB is the result of ‘combining’ the ith row of A with the jth or other adornments). 


column of B. 


Several ways have been developed to help apply the above definition. If 
you are good at mental arithmetic, one of the simplest ways is to lay your 
hands along the rows and columns of the matrices as you are multiplying 
them. For example, to calculate the top left element of the product matrix, 
place your hands as shown in Figure 2.1. (If you need to use pen and 
paper or a calculator, slips of paper can play the same role as the hands.) 


& 0.6 : | 
0.8 O04 O.7 0.4 


Figure 2.1 Laying hands on a matrix multiplication to outline rows and columns 


The corresponding elements can then be read off underneath the left hand 
and on the left of the right hand; in this case we obtain 


0.2 x 0.44 0.6 x 0.2 + 0.3 x 0.4. 


The next activity provides some practice in matrix multiplication. 
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Note that A? is the same size 
as A. 
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Activity 2.2 Calculating matrix products 


For each of the following pairs of matrices A and B, calculate the matrix 
product AB. 


For which of these pairs of matrices can you also find the matrix 
product BA? Calculate those which can be formed. For those which 
cannot be formed, say why this is so. 


wa-(42), B=() 34) 


Solutions are given on page 49. 


This way of combining matrices is called matrix multiplication because it 
shares many of the properties of multiplication of numbers. However, there 
is one important difference between the properties for matrices and for real 
numbers. The order in which you multiply two real numbers does not 
matter, since ab = ba for all real numbers a and b. This is not true for 
matrix multiplication. There are matrices A and B for which the product 
matrices AB and BA can be formed, but for which 


AB+BA. 


If you have not noticed this already, then look back at Activity 2.2(b), 
which shows that even if the product matrices are the same size, 
AB might not be equal to BA. 


Powers of matrices 
Consider the matrix 
2 1 
ea, 


Just as x? is defined as the product x x x, so we can define the square 
of A, written A?, as the matrix product AA; thus 


oa ate Te a Le 
BOA é Os '3 OF 1b 27° 
Not all matrices can be squared, however. The next activity considers the 


possibility of squaring a matrix. 


Activity 2.3 When does the square of a matrix exist? 


(a) Let M= (j ’ i) 


Why is it not possible to form the matrix product MM ? 
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(b) Use your answer to part (a) to formulate the condition on a matrix A 
under which it is possible to calculate its square A?. 


Solutions are given on page 49. 


Any matrix with the same number of rows as columns is called a 

; 21). 10 38). 
square matrix. So A = & a is square, but M = ¢ 9 is not. 
From the above activity, the condition for the square of a matrix A to 
exist is that A is a square matrix. 


We can also form the products (A?)A and A(A?), because A is a square 
matrix and A? is a matrix of the same size as A. For example, 


a a=(6 3)(3 o)=(a 6) 
aay=(3 0) (6 3) = (a 6): 


In this example (A?)A = A(A?), so we can write this product as AAA, 
without ambiguity. We denote AAA by A®; that is, the cube of A is 


A? = AAA. 


The cube of a matrix is always unambiguously defined because it does not 
matter which matrix multiplication is done first. In fact, more generally, it 
is true that 


(AB)C = A(BC) 


for any matrices A, B and C, in the sense that if the matrix product on 
either side of the equation exists, then so does the other and they are 
equal. One way of visualising this fact is to think about combining three 
networks. For example, consider the three networks shown in Figure 2.2, 
where network 1 is represented by the matrix C, network 2 is represented 
by the matrix B, and network 3 is represented by the matrix A. You can 
think of the entire network as network 1 and network 2 combined (BC), 
and then combining this with network 3 (A(BC)). Or you can think of it 
as network 2 and network 3 combined (AB), and then network 1 added on 
top ((AB)C). But changing the way you ‘think’ about the entire network 
does not change the physical network, so the flow through it is the same. 


2.2 Matrix addition 


Until now, all our efforts have been devoted to the motivation for, and 
mechanics of, multiplying matrices. This is because the importance of 
matrices arises from matrix multiplication. But matrices can also be 
added, and this is sometimes useful. (For example, vector addition has an 
important geometric interpretation.) 


Now matrix addition is introduced. Two matrices can be added only if 
they have the same size. To add two matrices, just add the corresponding 
elements. For example, 


Poof et B OY fite 228 48 
4 5 6 10 It 12/~ \4+10 5411 6412 
(8 10 12 
~\14 16 18)" 


Above, we found 


$17 2 
ar=(j 2). 


For any positive integer n, 
the nth power of a square 
matrix A, written A”, is 
defined in the same way, by 
repeated multiplication. 


This equation means that the 
matrix product ABC is 
unambiguous. 


network 


network 
3 y 


Figure 2.2. Combining three 
networks 


The geometric approach to 
vector addition is discussed in 
Chapter B3. 


21 


22 


CHAPTER B2 MODELLING WITH MATRICES 


The rule for adding matrices extends to the sum of any number of matrices 
of the same size: just add corresponding elements. In particular, this 
means that 


(A+B)+C=A+(B+O), 
so the expression A + B+ C is unambiguous. 


Now try matrix addition yourself. 


Activity 2.4 Adding matrices 


For each of the pairs of matrices below, calculate the matrix sums A + B 
and B+ A. 


wa-(3 a) a-(4 $) 
wae(28 9. we(b 2) 


40 5 6 20 
(c) A=/61 OO], B=[0 -1 4 
2 5 O01 


Solutions are given on page 49. 


Activity 2.4 demonstrates that order is unimportant in matrix addition; 
that is, A+ B= B+A for all matrices A and B of the same size. 


Matrix addition has an important special case, which is vector addition: 
the sum of two vectors with the same number of components is formed by 
adding corresponding components. For example, 


(sea) 


Now try adding vectors yourself. 


Activity 2.5 Adding vectors 


(a) Find the following vector sums. 


0 (88) (0) 
GEC) 
(°C) 


(b) Give an example of a pair of vectors which cannot be added together. 


Solutions are given on page 49. 
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Scalar multiplication of matrices 


Any matrix can be added to itself. Often A + A is written as 2A. This 
notation can be generalised to any real number & and any matrix A, so 
that kA is the matrix with each element equal to k times the 
corresponding element in A. For example, 


(3 jew Ge ‘) 
2 = $x0 $x(-1) 0 -3 


The matrix kA is called the scalar multiple of the matrix A by the real 
number k, and the operation of forming kA is called scalar In this context, the real 


multiplication of matrices. number k is often called a 


; F i scalar. 
Now try applying this operation for yourself. 


Activity 2.6 Scalar multiplication of matrices 


(a) In each case below, simplify the expression by performing the scalar 
multiplication. 


oa 2) ma?) wad 


(b) Simplify each of the following matrices by writing it as a scalar 
multiple of a matrix with integer entries. 


oP gs) @ () 


(iii) : ), where x is a real number. 


(c) Let A = € a) and B= € a Show that 


A(2B) = 2(AB). 


Solutions are given on page 50. 


Activity 2.6(c) illustrates the following useful result concerning scalar 
multiplication. If the matrix product AB exists, then 


A(kB) = (kA)B = k(AB), (2.2) 
where k is a real number. 


The operation of matrix subtraction is defined in terms of scalar 
multiplication and addition of matrices. The result of subtracting a 
matrix B from a matrix A, written A — B, is given by 
A—-B=A+4(-1)B, The matrix (—1)B is often 
written as —B, and is referred 
where (—1)B is the scalar multiple of B by —1. Thus, only matrices of the — to as the ‘negative of B’. 


same size can be subtracted, and the subtraction is performed by 
subtracting corresponding elements. The following example illustrates this 
operation: 


(<2 7)-(s -a)= (422s 22) = (00 *8): 
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A comma can be added if a 
subscript is ambiguous. For 
example, @1;2 could be 
written as a@11,2 OF @1,12, as 
appropriate. 


The following diagram might 


help you to remember this. 


col 1 col 2 
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Factorising 


When solving problems involving matrices, you may encounter a matrix 
expression of the form 


AB+AC. 
It can be shown that this expression can be factorised, to give 
AB+ AC = A(B+C). (2.3) 


This equation holds whenever the matrix products and sums on each side 
exist. 


Activity 2.7 asks you to verify equation (2.3) in a particular case, where 
B and C are vectors. 


Activity 2.7 Verifying a result 
1 -1 1 —1 
Let A= (5 1):B=(3) and C= ( |): Show that 
AB+AC=A(B+C). 


A solution is given on page 50. 


Element notation 


Sometimes it is useful to have a notation for the elements of a matrix. 
Consider the matrix 


($3) 


This is a 2 x 2 matrix, and its elements can be described in terms of row 
and column positions using subscripts: 


aii Ai2 
ie ( 
a21 29 


Note that the lower-case italic letter a has been used to denote the 
elements of the matrix A, to emphasise that they are numbers rather than 
matrices. The subscript is used to specify the position of the element, as 
follows: 


a, is the element in the first row and the first column (a1; = 1); 

aj is the element in the first row and the second column (a12 = 2); 
21 is the element in the second row and the first column (a2; = 3); 
22 is the element in the second row and the second column (dz2 = 4). 


The following activity tests your understanding of this notation. 
Activity 2.8 Interpreting element notation 
For the matrix 
7 2.5 0 
a ce 1 9) , 
write down the following elements. 


(a) bi4 (b) ba1 (c) bis 
Solutions are given on page 50. 
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As an example of how element notation can be used, here is a succinct 
statement of the rule for matrix addition. 


Matrix addition 


If C is the sum of two m x n matrices A and B, then C is an m x n 
matrix with elements 


Cig = iz + diy, 


where 1 = 1,2,...,m and 7 = 1,2,...,n. 


Element notation can also be used to state concisely the rule for matrix 
multiplication. To see how to do this, consider two 2 x 2 matrices 
A and B, and let AB = C. We can write 


AR = & a) ts ) = & | =i 
a21 G22 bo: dae C21 C22 
Now consider the top left element c,; of the product matrix. This element 


is formed by adding the products of the corresponding elements in the top 
row of A and the left-hand column of B: 


C11 = (a4, X by) + (iz X by). 


The same pattern of subscripts appears in the calculation of the other 
elements, and in general an element of the product matrix is given by 


Ci = (aia x bi;) + (ai2 x bo;), where i = WD j = 1,2. 


You can check this by multiplying out the above matrix product. This 
formula applies to any 2 x 2 matrix (or indeed any matrix product where 
the first matrix has two columns and the second matrix has two rows, for 
example, when multiplying a 3 x 2 matrix by a 2 x 3 matrix). 


For multiplying two 3 x 3 matrices, the rule is 
Cij = (ait x bi;) + (ai2 x bj) + (aiz x b3;), where i= 1, 2,3, J = 1, 2,3. 


This can be written, using the sigma notation for sums that was 
introduced in Chapter B1, as 


3 
Ci = ) Qik X Diy: 
k=1 


From this 3 x 3 case, we generalise to the case of multiplying any two 
matrices of appropriate sizes. 


Matrix multiplication 


If C is the product of an m x n matrix A and an n x p matrix B, 
then C is an m x p matrix with elements 


n 
Cij = S ik X dj, 
k=1 


where 7 = 1,2,...,m and j =1,2,...,p. 


See Chapter B1, 
Subsection 1.2. 


In 


Note that lower-case bold 


symbols are used for vectors. 


When writing vectors by 
hand, you should underline 
the symbol either with a 
straight line, v, or with a 
wavy line, v. 
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Element notation can also be used for the components of a vector. Since 
vectors always have only one column, the second subscript is always 1; for 
example, for a vector v with two components we have 


(m) 
v= : 
V21 


The above notation is correct, but cumbersome. We usually omit the 
second number, and write 


Summary of Section 2 


This section has introduced: 


© multiplication and addition of two suitable matrices A and B, written 
as AB and A +B, respectively; 


© scalar multiplication of a matrix A by a real number k, written as kA; 
© some properties of matrices: 

(AB)C = A(BO), 

A+B=B-+A, 

(A+B)+C=A+(B+O), 

A(kB) = (kA)B = k(AB), 

AB+AC = A(B+ OC); 


© the notation a;; for the element in the 7th row and jth column of the 
matrix A. 


Exercises for Section 2 


In Exercises 2.1 and 2.2, use the following matrices. 
I: <2 

2 0 4 6 5 1 4 2 
a=(7 7 B=(5 1 > cae 3 : p=(4 7 


Exercise 2.1 
(a) Which of the following matrix products can be formed? 


AB CA CD DA DB A’ D’ 


(b) Evaluate those matrix products in part (a) which can be formed. 


Exercise 2.2 


(a) Which of the following combinations of matrices can be formed? 
A+B B+C C+A C+B D+D A+A+A B-C 


(b) Evaluate those combinations in part (a) which can be formed. 
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Exercise 2.3 


Evaluate each of the following scalar multiplications. 
3 -l 1f2 95 3-5 
@5(9 3) 4(% 3) ©-3G 25) 


Exercise 2.4 


wl 


Evaluate each of the following. 


© (3)+(2) om (ot) - (08) 
© (3)+(3) © (4)+() 


Pa 


3 A population problem 


THE ECHO 


Figure 3.1 A newspaper 
headline 


FF speci 


Data adapted from 
Demographic Yearbook 1991, 
© United Nations 1992. 
This book uses 15 as the age 
group boundary for 
adulthood, and 65 as the age 
of retirement, so we adopt 
these conventions too. 
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In Chapter B1 you studied the long-term behaviour of populations, using 
recurrence sequences to explore the ways in which a population increased 
or decreased — or perhaps remained the same. But many of the questions 
asked about populations are concerned not so much with the total 
population as with its structure. For example, underlying the headline 
question in Figure 3.1 is a concern about the proportions of pensionable 
and working people in the population, and how these groups may vary in 
the future. This suggests that a model to address this question will need to 
divide the total population into a number of subpopulations, and then to 
explore the ways in which these subpopulations behave. That is the 
approach adopted in this chapter. 


3.1 Population modelling 


In any human population, there are two major groups which are not part 
of the workforce: the juveniles who have not yet entered the workforce and 
the elderly who have retired. So to begin to investigate the question of 
whether the potential workforce is big enough to support the whole 
population, we shall divide it into three broad age groups — juveniles, 
workers and elderly — and investigate the proportions of each of these in 
the total population. Table 3.1 shows data for the human population of 
the UK in June 1990. (All values in the table are correct to the number of 
decimal places given.) 


Table 3.1 UK population in June 1990: juveniles, workers and elderly 


Juveniles Workers Elderly 
(aged < 15) (aged 15-64) (aged > 65) 
Subpopulation (millions) 10.92 37.50 8.99 
Birth rate 0.0000 0.0213 0.0000 
Death rate 0.0008 0.0031 0.0575 


The birth rates quoted in Table 3.1 are the annual proportionate birth 
rates, as defined in Chapter B1. For example, the number 0.0213, which 
we shall call the ‘worker birth rate’, is calculated as follows: 


number of babies born with worker parents 


worker birth rate = (3.1) 


total number of workers 
The death rates quoted in Table 3.1 are the annual proportionate death 
rates, and are calculated similarly. Test your understanding of these data 
by trying the following activity. 


Activity 3.1 Understanding the data 


(a) Use the data in Table 3.1 to estimate the number of babies born 
between June 1990 and June 1991. 


(b) Use the data in Table 3.1 to estimate whether more workers or elderly 
died in the period from June 1990 to June 1991. 


Solutions are given on page 50. 


SECTION 3 A POPULATION PROBLEM 


We have divided the population into three age groups in order to start to 
address the question posed in Figure 3.1. However, to establish a technique 
for modelling such subdivided populations, we look first at the simpler 
case where there are just two subpopulations — juveniles and adults. In 
Section 4 you will analyse the original question on the computer. The 
simpler question that we shall investigate in this section is the following. 


Will there be a constant proportion of juveniles entering the potential 
workforce, or will we have an ever-increasing proportion of adults? 


The data appropriate for this simplified problem are shown in Table 3.2. 


Table 3.2. UK population in June 1990: juveniles and adults 


Juveniles Adults 

(aged < 15) (aged > 15) 
Subpopulation (millions) 10.92 46.49 
Birth rate 0.0000 0.0172 
Death rate 0.0008 0.0136 


Note that the data in Table 3.2 are calculated directly from the data in 
Table 3.1. So, for example, 


janie teats total number of births to workers and to elderly 


total number of workers and elderly 
__ 0.0213 x 37.50 + 0.0000 x 8.99 
7 37.50 + 8.99 


~ 0.0172. 


The death rate for adults is calculated in a similar way, from the death 
rates for workers and for elderly. 


Now we are going to create a mathematical model to estimate the numbers 
of juveniles and adults in the years following 1990. The first task is to 
assign variables to the quantities in the model. 


Table 3.2 gives figures for June 1990, and we shall take this date as the 
start of year 0, with year n starting n years after that date. Suppose that 
at the start of year n, the number of juveniles is J;,, and the number of 
adults is A, (in millions). In terms of these variables, Ag is the number of 
millions of adults in June 1990 as given in Table 3.2, so Ag = 46.49. 
Similarly, Table 3.2 gives Jo = 10.92. 


How do we calculate J; and A,? In Chapter B1 the (single) population 
estimate for year 1 depended on only the (single) population estimate for 
year 0. Here, because a juvenile either becomes an adult or stays a 
juvenile, and because an adult stays an adult but can produce new 
juveniles (babies), we can expect that each of the two subpopulation 
estimates for year 1 will depend on both of the subpopulation estimates for 
year 0. This interdependency can be represented by the network shown in 
Figure 3.2. There are four pipes shown in Figure 3.2, which represent the 
‘flow’ from one subpopulation to another. The ‘quantity of flow’ down each 
pipe can be given an interpretation in this current context, which is 
summarised in Table 3.3. Using this interpretation, we can determine the 
labels for the pipes. 


Jo Ao 
y y 
Jy Ai 


Figure 3.2. A network which 
shows the interrelationships 
between subpopulations in 
successive years 


7118] 
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Table 3.3 Interpretation of the pipes in Figure 3.2 


From ‘To _ Interpretation 


Jo J, Proportion of juveniles who neither die nor reach age 15 in year 0 
Jo A, Proportion of surviving juveniles who reach age 15 in year 0 

Ao J, Proportionate birth rate for adults in year 0 

Ap A, Proportion of adults who do not die in year 0 


(Note that it is convenient to use Jy, J;, Ag and A; as names for the nodes 
in the network as well as for estimates of subpopulation sizes.) 


Consider the pipe from Jp to J,;. There are Jy juveniles at the start of 
year 0. From the data given in Table 3.2, a proportion 0.0008 of these die 
during the year, so a proportion 1 — 0.0008 = 0.9992 are alive at the end of 
the year. Some of this number will have reached the age of 15 during the 
year and become adults. We shall assume that iz of the surviving juveniles 
become adults. (This is the case if there is a uniform spread of ages within 
the group.) So the proportion of juveniles still alive and still juveniles at 
the end of the year is iz x 0.9992 = 0.9326 to four decimal places. Thus 
the number that labels the pipe from Jp to J; is 0.9326. 

Now try to calculate the numbers that label the other pipes in the 
network. 


Activity 3.2 Labelling the pipes 


There are four pipes in the network shown in Figure 3.2, and each has an 
interpretation given in Table 3.3. The number labelling the pipe from Jp 
to J, was calculated above, using the assumption that one fifteenth of 
surviving juveniles become adults each year. Now calculate the labels for 
the other pipes (correct to four decimal places). 


(a) Calculate the number that labels the pipe from Ag to J;. 

(b) Calculate the number that labels the pipe from Ag to Aj. 

(c) Calculate the number that labels the pipe from Jy to Aj. 
(Hint: Use the ‘’ assumption.) 


Solutions are given on page 50. 


The network model of the changes in the two subpopulations from year 0 

to year 1 can be converted into a matrix model using the correspondence 

established in Subsection 1.2. So the estimates for the two subpopulations 
in year 1 can be calculated by 


J,\ — (0.9326 0.0172 Jo 

A,] ~~ \0.0666 0.9864] \ Ap /° 
If we assume in our model that the birth and death rates are constant 
(that is, the same each year), then the same matrix can be used to 


calculate the subpopulation estimates for subsequent years. Thus we have 
the matrix model 


Ing \ — (0.9326 0.0172 In _ 
a) - oie ee) (3°) ess) 34) 
Now let us pause and consider some simplifying assumptions that were 
made or could have been made in arriving at these equations. 
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Activity 3.3 Assumptions in the model 


For each assumption in the list below, state whether it was needed to 
create the model in equation (3.2), and comment on whether you think the 
assumption is reasonable. 


e The birth and death rates are constant. 
e No juveniles give birth. 

e One fifteenth of surviving juveniles become adults each year. Fy, An 
e There is no emigration or immigration. 


Solutions are given on page 51. 


0.9326 Y¥ y 0.9864 


Since the birth and death rates are assumed constant in this model, the 
network diagram in Figure 3.2 can be generalised to year n and year n+ 1, 


as shown in Figure 3.3. This figure includes the calculated labels. Fath An+1 
Now you are asked to use the matrix model to calculate the population Figure 3.3 
estimates for year 1 and year 2. Generalising Figure 3.2 


Activity 3.4 Calculating 


Use equation (3.2) and the data in Table 3.2 to estimate the values of Fo) 
Ji, Al, J and Ao. — 4 


State what each of these numbers represents. 


Solutions are given on page 51. 


By repeating the type of calculation in Activity 3.4, we obtain predictions 
for the seven years following 1990, shown in Table 3.4. (In this table and 
related ones, each prediction is for June of the given year, and n = 0 
corresponds to 1990.) 


Table 3.4 Predictions from the model (populations in millions) There is no need for you to 
check the predictions for 
1993-1997. 


1990 1991 1992 1993 1994 1995 1996 1997 


Juveniles, J, 10.92 10.98 11.04 11.10 11.16 11.21 11.27 11.32 
Adults, Ay, 46.49 46.59 46.68 46.78 46.89 46.99 47.10 47.21 


Looking at the results shown in Table 3.4, it can be seen that the numbers F iiterpretily 
for both subpopulations are increasing. First we investigate the increase in ’ lee 
the total population, by defining 7}, to be the total population at the 

beginning of year n, that is, T;, = J, + An. Now the variation in the total 

population can be seen by calculating T;,, as shown in the first row of 

Table 3.5. 


Table 3.5 Predictions of total population from the model 


1990 1991 1992 1993 1994 1995 1996 1997 


Total, T, (millions) 57.41 57.57 57.73 57.89 58.05 58.21 58.37 58.53 
Ratio T,,/T,-1 1.003 1.003 1.003 1.003 1.003 1.003 1.003 


From Table 3.5 it can be seen clearly that, according to the model, the 
population is increasing geometrically, with the ratio between successive 
years equal to 1.003 (to three decimal places). This means that the annual 


proportionate growth rate 
is 0.003. 
Sik 


Note that the two proportions 


add to 1, because 


Jn 
ee 


Tn 
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Apart from this general trend upwards, is there any other information to 
be gleaned from Table 3.4? To answer this, it is best to eliminate the 
general upward trend in each subpopulation, by considering the proportion 
of each in the total population. So we divide both subpopulation figures by 
the total population for the corresponding year. The results are tabulated 
in Table 3.6. 


Table 3.6 Predictions of the subpopulation proportions from the model 


1990 1991 1992 1993 1994 1995 1996 1997 
Proportion of 0.190 0.191 0.191 0.192 0.192 0.193 0.193 0.193 
juveniles, J;,/Ty, 
Proportion of 0.810 0.809 0.809 0.808 0.808 0.807 0.807 0.807 
adults, An/Tp 


Looking at Table 3.6, you can see that the proportion of juveniles is 
gradually increasing. In Section 4, by plotting a graph on a computer, you 
will see that the rate of increase is slowing down. Also, by using a 
computer to calculate the predictions for the first fifty years, you will show 
that the proportion of juveniles approaches the limit 0.197, so the model 
predicts that eventually 19.7% of the population will be juveniles. The 
proportion of adults approaches the limit 0.803, so the model predicts that 
80.3% of the population will be adults. So the answer to the question 
posed on page 29 is that the model predicts a stable situation where a 
constant proportion of juveniles enters the workforce each year, namely 
19.7/15 = 1.3 per cent of the population. 


So the proportions of juveniles and adults eventually settle down into 
constant proportions. You have already seen something analogous to this 
at the end of the video sequence. The last experiment involved cycling 
water repeatedly through a two-input, two-output network. It was found 
that the proportions quickly settled down. 


To evaluate a model, more data are needed. The Demographic Yearbook 
(1997 edition) gives 11.36 million juveniles and 47.44 million adults 

in 1996. Comparing these figures with the values of Jg = 11.27 and 

Ag = 47.10 given in Table 3.4, the agreement is pretty good but not 
perfect, which suggests that the model should be refined. You will return 
to the more detailed three-subpopulation model in the next section, and 
will use it to answer the question ‘Who will support us in our old age?’ 
that was posed at the beginning of this section. 


3.2 Closed-form solution 


The matrix model of the subpopulations given by equation (3.2) has a 
closed-form solution. We derive this here. 


Equation (3.2) can be expressed more succinctly as 
Pr4i — Mp,, (n — 0, 1s 2, ee as (fa) 


; _ {0.9326 0.0172 : dy 
where M is the matrix ee one) and p,, is the vector Ce 


Now consider how successive estimates of the subpopulations are obtained. 
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The population estimate for year 1 is obtained directly from the initial 
data, as p, = Mp,. The population estimate for year 2 is obtained from 
the population estimate for year 1, as p, = Mp,. Substituting for p, from 
the previous equation gives the population estimate for year 2 in terms of 
the initial data: p, = M(Mp,) = M’p,. Similarly, the estimate for year 3 
can be obtained from the initial data as p; = Mp, = M(M?’p,) = M°pp. 
Continuing in this way, we obtain the estimate for the population in year n 
in terms of the initial data, as 


Pr = M"po. (3.4) 


This is the closed-form solution for the population problem. 


As it stands, the closed-form solution (3.4) is no improvement on the 
recurrence system (3.3) as in fact it takes more arithmetic operations to 
compute the former than the latter. The closed-form solution is useful, 
however, because there is a quick way of calculating powers of a matrix. 
This quick method is beyond the scope of this course, but is covered in 
second-level courses. 


Summary of Section 3 


This section has introduced a matrix model for exploring the ways in 
which two interdependent subpopulations behave, in order to gain some 
insight into the long-term behaviour of the structure of a population. 


Exercise for Section 3 


Exercise 3.1 


The aim of this question is to develop a model of the population of a 
fictitious developing country. Table 3.7 gives the birth and death rates for 
two subpopulations of that country. 


Table 3.7 Population of a fictitious developing country 
in January 2000: juveniles and adults 


Juveniles Adults 

(aged < 15) (aged > 15) 
Subpopulation (millions) 47.85 614.37 
Birth rate 0.0000 0.0520 
Death rate 0.0139 0.0231 


(a) Draw a network diagram similar to Figure 3.2, and label the pipes 
using the interpretations given in Table 3.3 and the data given in 
Table 3.7. 


(b) Write down the corresponding matrix model for the population. 


You saw this approach for 
obtaining a closed-form 
solution in Chapter Al. 
There you saw that the 
recurrence equation 


In4+1 = k&n, 


where x, and k are real 
numbers, has the closed form 


tn = kh” xo. 
This was put to use in 
Chapter B1, in obtaining a 
closed-form solution for the 


exponential model of 
population change. 
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4 Computing with matrices 


cs 


Note that there are no pipes 
from J, to E,+1 and 

E,, to Wy41: juveniles cannot 
become elderly in one year, 
and the elderly cannot 
become younger, to become 
workers. There is also no link 
from E,, to Jn41, on the 
assumption of a zero birth 
rate for the elderly. 
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In this section you will be using Mathcad to create, add and multiply 
matrices. This work will increase your familiarity with the way matrices 
behave; also, it will help you to distinguish between circumstances when it 
is worth using Mathcad to perform matrix calculations, and when it would 
be more efficient to carry out the calculations on paper. 


You will also investigate further the two-subpopulation model studied in 
Section 3, and solve the three-subpopulation model set up there. This 
model will enable you to address the question of whether in the long term 
there will be enough workers to support the ageing population. 


The three-subpopulation model is encapsulated in the following network 
diagram and table, where at the start of year n, J, is the number of 
juveniles, W,, is the number of workers, and E,, is the number of elderly. 


0.9326 Y VY 0.9425 


Wh +1 


Fyt1 


Figure 4.1 The interrelationships between the subpopulations in successive years 


Table 4.1 Interpretation of the pipes in Figure 4.1 


From To Interpretation (for year 7) 


J Jn+1 Proportion of surviving juveniles who do not reach 15 
Wr  Wnr+i Proportion of surviving workers who do not reach 65 
FE, Ens Proportion of elderly who do not die 

Wn Jn+1  Proportionate birth rate for workers 

Ji Wn+1 Proportion of surviving juveniles who reach 15 

Wr, En+i Proportion of surviving workers who reach 65 


The pipe labels in Figure 4.1 are calculated by using the interpretations in 
Table 4.1, as in Section 3 for the two-subpopulation model. 


Refer to Computer Book B for the work in this section. 


Summary of Section 4 


This section provided practice in using Mathcad to create, add and 
multiply matrices. The population problems described in the previous 
section were studied in more detail. 


5 Simultaneous linear equations and matrices 


In this section we explore the close connection between simultaneous linear 
equations and matrices. In particular, you will see that simultaneous linear 
equations can be solved by using only matrix methods. We start the 
section by exploring a little more of the theory of matrices, and then apply 
it to solve simultaneous linear equations. In order to keep the arithmetic 
to a minimum, this section will be concerned mostly with 2 x 2 matrices. 


5.1 Identity matrices 


An identity matrix is a special type of matrix: multiplying a matrix by an 
identity matrix leaves the matrix unchanged, if the product exists. Do the 
following activity to see how this works. 


Activity 5.1 A case of identity 


3.4 2 7 -3 1 0 
Let A=(_5 ine; l ) and 1= (4 a) 


Evaluate each of the following matrix products, where possible: 
AI, IA, BI, IB. 


Solutions are given on page 51. 


Let I be the matrix € a and let A be any matrix for which the 


product AI can be formed. (Thus A must have 2 columns.) Then it is 
always true that 


AI=A. 
What is more, if B is any matrix for which the product IB can be formed 
(that is, B is any matrix with 2 rows), then 

IB = B. 
This matrix I is called the 2 x 2 identity matrix. 


Note that if A is a 2 x 2 matrix, then the products AI and IA can both 
be formed since both conditions above concerning rows and columns are 
met. So, for a 2 x 2 matrix A, we have 


AI=IA=A. 
There is only one 2 x 2 identity matrix (the matrix I used above), but 
there are also identity matrices of other sizes; for example, The identity matrices of sizes 
2 x 2 and 3 x 3 are sometimes 
1 0 0 ; 
1-101 0 written as I, and Is, 
~~ respectively. 
001 


is the 3 x 3 identity matrix. This also has the properties AI = A and 
IB = B, for suitable matrices A and B (see Activity 5.2 for a verification A must have 3 columns, and 
of this assertion). B must have 3 rows. 


SD 


This situation is analogous to 
the arithmetic of numbers. 
For example, dividing by 2 is 
equivalent to multiplying by 
1 


5 = 271, which is the inverse 


of 2. 
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Each identity matrix is a square matrix with ones down the leading 
(or main) diagonal (that is, the diagonal of the matrix which starts at the 
top left and ends at the bottom right) and every other element zero. 


Activity 5.2 Identity matrices 


Suppose that I is the 3 x 3 identity matrix and b is a vector with 
3 components. 


(a) According to the assertion in the text above, what is the matrix 
product Ib? 


(b) Now write b as 


and calculate Ib, thus verifying that the assertion is correct in this 
case. 


Solutions are given on page 51. 


5.2 Inverse of a matrix 


You may have realised by now that we have not divided a matrix by a 
matrix — there is no useful way of doing this. However, in some 
circumstances, multiplying by an inverse matrix suffices in situations 
where we would want to divide. This subsection describes what we mean 
by an inverse matrix and how to calculate the inverse of a 2 x 2 matrix. 
We start by looking at some special pairs of matrices. 


Activity 5.3 Linked pairs 


(a) Multiply the following pairs of matrices together. In each case, try 
reversing the order of multiplication (BA instead of AB, for example), 
and see whether this affects the outcome. 


w a=(5 2) B=(3 8) 
a) C=(7 3), D=(4 9) 


(iii) P= wy a=({ 3) 


(b) Describe the outcome of your calculations in part (a). 


Solutions are given on page 51. 


Suppose that A and B are any two 2 x 2 matrices with the property that 
AB =I=BA, 
where I is the 2 x 2 identity matrix (such as the matrices A and B in 


Activity 5.3). Then we say that B is the inverse of A. We write the 
inverse of A as A~', so we have 


AA =T=A“*A, 
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This equation shows that if A~‘ is the inverse of A, then A must be the 
inverse of A~!. For the particular matrices in Activity 5.3, you showed 
that A = B~! and B= A“!, and that C= D1 and D=C"!. The 
matrices P and Q are not the inverses of each other, however, since their 
product is not an identity matrix. 


Even though the matrices P and Q are not inverses of each other, they are 
special in that the product of P and Q is a scalar multiple of the identity 
matrix, namely 31. Since PQ = 3I, 


1(PQ) = 4(31) =L 
Using the result that A(KB) = k(AB), given in equation (2.2), we have 
P(;Q) = ;(PQ) =I. 
Similarly, ({Q)P =I. So the inverse of P is obtained by scalar multiplying 
Q by re This ‘trick’ of scalar multiplying matrices appropriately can be 


used to find a general formula for the inverse of a 2 x 2 matrix. The next 
activity prepares the ground. 


Activity 5.4 Deriving a formula for the inverse of a 2 x 2 
matrix 


A general 2 x 2 matrix can be written in the form é a where a, b, ¢ 
and d are real numbers. Find the product 
6 4 ( d a 
c d —c aj- 
A solution is given on page 52. 


Comment 


Element notation has not been used in this example, since it would be 
cumbersome to do so. 


The product 


(oa) (ea) 


in Activity 5.4 yields a scalar multiple of the identity matrix, namely 
(ad — bc)I, so scalar multiplying the second matrix above by 1/(ad — bc) 
gives a formula for the inverse of a general 2 x 2 matrix. 


Inverse of a 2 X 2 matrix 


The inverse of the general 2 x 2 matrix (< )) is given by 


ad — bc & _ ) provided that ad — bc of 0. 


PQ = QP =31 


Similarly, 3P is the inverse 


of Q. 
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The following example shows how this formula can be used to compute the 
inverse of a matrix. 


Example 5.1 Finding an inverse 


Find the inverse of the matrix 
4 7 
a=(3 6) 
and check that your answer is the inverse. 


Solution 


By comparing the matrix A given in the question with the general 2 x 2 
matrix, we see that a= 4, b=7,c=3 andd=6. So 


ad — bc = (4 x 6) — (7 x 8) = 24-21 =3, 


and the inverse of A is 


Att=l 6 =f = 2 = 
3\-3 A a 2 


To check that this 7s the required inverse matrix, we multiply it by A: 
( 2 = ¢ a 2 wee 2x 7+ (-2) S) 
-1 ¢/)\3 6 (-1)x4+$x3 (-1)x7+$x6 
_ fi O 
SO aL" 


This product is the identity matrix, so we have checked that A~! is the 
inverse of A. 


It is useful to be able to refer to the quantity ad — bc: it is called the 
Square matrices of size bigger determinant of the 2 x 2 matrix 


than 2 x 2 also have a 
determinants, but they are A= ( ) ; 


not discussed in this course. cd 
and is often written det A. Thus 
det A = ad — be. 


If the determinant of a matrix is zero, then the above formula for the 
inverse of the matrix fails (because division by zero is not defined). Not 
only does the formula fail but, in fact, a matrix with zero determinant does 
not have an inverse, as you will see geometrically in the next subsection. 
This role of the determinant is summarised in the following test. 


Determinant Test 


Other terms in widespread If the determinant of a matrix A is not zero, then A has an inverse 
use are non-singular instead and we say that A is invertible. 
of invertible and singular 


. : : If the determinant of a matrix A is zero, then A does not have an 
instead of non-invertible. 


inverse and we say that A is non-invertible. 
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3.4 ; 
6 5) has determinant 


det B=3x8-4x6=24-24=0. 


For example, the matrix B = ( 


Thus B is non-invertible; that is, B~' does not exist. 


The next activity offers practice in identifying 2 x 2 matrices that have an 
inverse and in finding inverse matrices. 


Activity 5.5 Finding inverses 


Use the Determinant Test to classify each of the following matrices as 
invertible or non-invertible. For each of the invertible matrices, write down 
its inverse. 


@A=(2 5) © B=(5 5) 


@e=(4 3] (@) D=(95 is) 


Solutions are given on page 52. 


5.3 Simultaneous linear equations 


In Chapter AO, two methods (substitution and elimination) for solving 
pairs of simultaneous linear equations were described. Here a third 
method, which uses matrices, is described. 


An example of a pair of simultaneous linear equations in the two variables 
x and y is 


0.42 + 0.2y = 4, 


5.1 
0.6x + 0.8y = 11. eH) 


Each of these equations represents a straight line in the (z, y)-plane. You 

may not recognise the above equations as representing straight lines, but 

they can both be rearranged into the more familiar form y= ma-+c. This The representation of a 

can be done by subtracting the x-term from both sides and then dividing straight line by an equation of 


by the coefficient of y. So the first equation above becomes the form y = mx + c was 
discussed in Chapter A2, 
4 — 0.4% . 
y= 2 = —2x + 20, Subsection 1.1. 
and the second equation becomes 
11 — 0.6z 3 55 
u=—"08 aot 


These two straight lines are shown in Figure 5.1, overleaf. 
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vA 


Figure 5.1 The two straight lines which correspond to equations (5.1) 


Equations (5.1) can be thought of as one equation connecting two vectors, 
rather than two equations connecting numbers, by writing each equation 
as one component of a vector, as follows: 


0.4¢+0.2y\ (4 

0.62+0.8y/ \11/)° 
Re-writing the left-hand side of this equation using a matrix—vector 
product gives the matrix form of equations (5.1): 


Ce elta: (5.2) 


We shall need to refer frequently to the 2 x 2 matrix that appears in the 
above equation, so we call it the coefficient matrix of the simultaneous 
linear equations (5.1). Note that the elements of the coefficient matrix are 
the same numbers arranged in the same order as the coefficients on the 
left-hand sides of equations (5.1). This will always be the case for the 
coefficient matrix of any pair of simultaneous linear equations in which the 
variables appearing on the left-hand sides are written down in the same 
order in each equation. 


The following activity is designed to give you practice in writing down the 
matrix form of a pair of simultaneous linear equations. 


Activity 5.6 Matrix form of simultaneous linear equations 


Rewrite each of the following pairs of equations in matrix form. 


(a) 22+ 3y = 38 (b) 2a =6 (c) 2a — 4y = 18 
x+4y=-—-1 3x + 6y = 15 49+ 3y=-1 


Solutions are given on page 52. 
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We now return to equation (5.2), the matrix form of equations (5.1), and 
proceed to solve the equations. Let A be the coefficient matrix, that is, 


0.4 0.2 
nS ty i) ; 
The determinant of A is det A = 0.2, so the inverse of A is 
het + 0.8 —0.2 _ 4 —-l 
0.2\-0.6 04 —-3 2/)° 


Now multiply both sides of equation (5.2) by A~', to obtain 


Ak (5) = Aq (1) (5.3) 


By definition, A~'A =I, since A~! is the inverse of A, and 
I (*) = (5). since multiplying a matrix by an identity matrix leaves 
that matrix unchanged. Hence, on simplifying the left-hand side and 


evaluating the right-hand side of equation (5.3), we obtain 
G\ 4 —-1 4\ (5 
yy Le OP yp Vaio? 

In other words, 


x=5 and y=10. 


The example above illustrates a general method for solving simultaneous 
linear equations which can be summarised as follows. 


Solving a pair of simultaneous linear equations using 
matrices 


First write the simultaneous linear equations in matrix form Ax = b, 
where A is the coefficient matrix, x is the corresponding vector of 
variables, and b is the vector with components equal to the 
corresponding right-hand sides of the equations. 


If the matrix A is invertible, then the solution is given by 


x= A'b. 


The next activity provides practice in using this method to solve 
simultaneous linear equations. 


Activity 5.7 Solving simultaneous linear equations 


Use matrices to solve the following pairs of simultaneous linear equations. 


(Hint: You have already written each pair of equations in matrix form in 
Activity 5.6.) 


(a) 24+ 3y=3 (b) 26=6 (c) 2a—2y=18 
e+4y=-1 3x2 + 6y = 15 Ay + dy =| 


Solutions are given on page 52. 


det A = 0.4 x 0.8 — 0.2 x 0.6 
= 0.32 — 0.12 = 0.2. 


This solution can be checked 
by substituting the values for 
x and y into equations (5.1). 


41 


An equation in n variables is 
linear if the largest power of 
any variable is 1, and no 
products of variables are 
present. For example, 


e+2yt+z=83 
is linear, but neither 
r+y’=2 
nor 
a+tyz=1 
is linear. 


42 


CHAPTER B2 MODELLING WITH MATRICES 


Each equation of a pair of simultaneous linear equations represents a 
straight line. So the solution to a pair of simultaneous linear equations 
must correspond to a point that is simultaneously on both lines, that is, 
the point of intersection. This is illustrated in Figure 5.2, which shows the 
lines which represent the pairs of equations that you solved in Activity 5.7. 


YA YA 
3a + 6y = 15 
> > 
rc 1 zr 
22 = 6 
(a) (b) 
UA 


RV 


(c) 
Figure 5.2. Points of intersection 


The matrix method of solution can be used for larger systems of linear 
equations. In general, a system of nm simultaneous linear equations in 

n unknowns can be written as a single equation in matrix form in which a 
known n-component vector b is equal to an unknown vector x multiplied 
by an n x n coefficient matrix A: 


Ax =b. 


If A is invertible, then this equation may be solved by multiplying both 
sides by the inverse of A, giving x = A~'b. The drawback of this method 
is that finding an inverse is a lengthy calculation for any matrix larger 
than the 2 x 2 matrices considered here. 


Exceptional cases 


Most pairs of simultaneous linear equations are like the ones considered 
above and have a unique solution. However, this is not always the case, 
and it is the exceptions that concern us next. 


We know that some 2 x 2 matrices do not have inverses. In the next 
activity you are asked to consider what this means in the context of 
simultaneous linear equations. 
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Activity 5.8 Non-invertible coefficient matrices 


Rewrite the following pairs of equations in matrix form. Show that in each 
case the coefficient matrix does not have an inverse. Try to solve the 
equations using the elimination method: what happens? 
(a) 2¢e-—3y=5 (b) a«+2y=-6 

—4x + 6y =7 —3x —6y = 18 


Solutions are given on page 53. 


You saw in the previous activity two pairs of equations for which the 
coefficient matrix did not have an inverse, so the matrix method of finding 
solutions would not work. You also saw that the elimination method 
provided no solutions. In fact, no method for finding a unique solution will 
work, because neither pair of equations has a unique solution. To 
demonstrate this visually, the straight lines that correspond to the 
equations have been plotted in Figure 5.3. 


—4ar 4+ 6y = 7 


av 
av 


22 — 3y =5 a+2y=—6 
—3u — 6y = 18 


(a) (b) 
Figure 5.3 Exceptional cases 


As you can see from Figure 5.3, there are two cases for which there is no 
unique solution. In Figure 5.3(a) the two lines are parallel, so they have no 
point of intersection. This shows that the equations have no solution. 

In Figure 5.3(b) both equations describe the same straight line, namely 

y = —x/2— 3. So the ‘two lines’ intersect at every point along the line, 
and there are infinitely many solutions. 


If the coefficient matrix is non-invertible, then one of the two possibilities 
described above occurs: either there is no solution or there are infinitely 
many solutions. In this course we do not elaborate on the methods of 
distinguishing between these two possibilities, beyond the graphical 
method used above. 


Final remarks 


The discussion in this subsection has shown how to solve pairs of 
simultaneous linear equations by finding the inverse of the coefficient 
matrix. The methods described in Chapter AO are preferred for hand 
calculations, because these methods involve fewer arithmetic operations 
than the matrix method. Note, however, that you may be specifically asked 
to use the matrix method for hand calculations in the assessment of this 
course, because it provides a good test of your understanding of matrices. 


A method for distinguishing 
between the two possibilities 
is described in second-level 
mathematics courses. The 
method is based on the 
‘elimination’ method 
described in Chapter AO, and 
is called Gaussian elimination. 
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For computer calculations, a method similar to the ‘elimination’ hand 
calculation method (described in Chapter AO) is quicker than finding the 
inverse of the coefficient matrix. 


So for what is the matrix form of the solution useful? The main use is in 
developing the theory of simultaneous linear equations. As the first payoff 
from this theory, we can say that simultaneous linear equations have a 
unique solution if and only if the coefficient matrix is invertible. Using the 
Determinant Test for invertibility then gives a quick test for the existence 
of a unique solution (or, equivalently, for a pair of lines having a unique 
point of intersection). Try this now in the following activity. 


Activity 5.9 Determinant Test for a unique solution 


Which of the following pairs of simultaneous linear equations has a unique 
solution? (You are not asked to find the solution.) 


(a) 5¢—3y=5 (b) 38a2-— y=-6 (c) 6” —2y=1 
—x+4y=7 —9r7+ 3y = 18 3x2 — Ty =—-1 


Solutions are given on page 53. 


Now would be a good time to listen to the optional band of CDA5505, 
which describes some applications of matrices. 


Listen to CDA5505 (Track 8), band 3, ‘Applications of matrices’. 


Summary of Section 5 


This section has introduced: 

© identity matrices; 

© the Determinant Test for identifying invertible matrices; 
© the formula for the inverse of an invertible 2 x 2 matrix; 
© 


expressing a pair of simultaneous linear equations in matrix form as 


Ax =b; 


© amethod for solving a pair of simultaneous linear equations Ax = b, 
in the case where A is invertible. 


Exercises for Section 5 


Exercise 5.1 


Find the inverse of each of the following matrices. 


@) A=(7 4) (») B= (5 ) ©) e=( 


Exercise 5.2 


Solve the following pairs of simultaneous linear equations using matrix 

methods. 

(a) —32 —5y=1 (b) 64+ 2y=4 (c) 4a + 2y = 2 
z+4y=2 32 — y=10 z+ 3y = —3 


Summary of Chapter B2 


In this chapter you met the concept of matrices. You have used matrices 
to investigate a population case-study and also to solve pairs of 
simultaneous linear equations. 


Learning outcomes 


You have been working towards the following learning outcomes. 


Terms to know and use 


Network, vector, component of a vector, matrix, element of a matrix, 
row of a matrix, column of a matrix, size of a matrix, matrix addition, 
matrix multiplication, square matrix, powers of a matrix, 

scalar multiplication of matrices, identity matrix, inverse of a matrix, 
invertible matrix, non-invertible matrix, determinant of a matrix. 


Symbols and notation to know and use 
pe eo) A= & ), Av, AB, A", A+B, 
V2 G21 422 


kA, A—B,I, Aq, detA 


Mathematical skills 
© Represent networks by matrices. 
© Describe a matrix in terms of the numbers of its rows and columns. 


© Write down an example of an m x n matrix for given values of 
m and n. 


Evaluate matrix products. 
Evaluate matrix sums. 
Evaluate the square of a matrix. 


Evaluate the result of scalar multiplying a matrix. 


oOo 9 9 2 


Use element notation to describe matrix elements in terms of their row 
and column positions. 


© Write down conditions under which it is possible to form matrix 
products, sums and powers. 


© Decide whether a 2 x 2 matrix has an inverse, and evaluate the inverse 
(if it exists). 


© Write a pair of simultaneous linear equations in matrix form, and solve 
them using matrix methods where possible. 


Mathcad skills 

© Create matrices. 

© Add and multiply matrices. 

© Scalar multiply matrices and evaluate powers of matrices. 
© 


Calculate a sequence of vectors using matrix—vector multiplication. 
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Solutions to Activities 


Solution 1.1 


(a) 


In Figure 1.1 the label on the pipe from node B 
to node U is 0.3. So 0.3 litres of water are 
output at node U. 


Similarly, by reading the labels on the other 
pipes from node B, we have 0.1 litres of water 
output at node V, and 0.6 litres of water output 
at node W. 


First consider inputting 1 litre of water at 

node A. Then, as described in the text 
preceding the activity, 0.4 litres of water will be 
output at node U. Now add another 1 litre of 
water to node A. This will add another 0.4 litres 
of water to the 0.4 litres already output at 

node U. So, all together, 2 x 0.4 = 0.8 litres of 
water will be output at node U. 


Similarly, the outputs at the other two nodes 
will double. So 2 x 0.2 = 0.4 litres of water will 
be output at node V, and 2 x 0.4 = 0.8 litres of 
water will be output at node W. 


Solution 1.2 


(a) 


Adding together the contribution from each 
input gives: 

0.4 + (2 x 0.3) = 1 litre of water output at 
node U; 

0.2 + (2 x 0.1) = 0.4 litres of water output at 
node V; 

0.4 + (2 x 0.6) = 1.6 litres of water output at 
node W. 


First consider the amount of water reaching 
node U from node A. 1 litre input gives 0.4 litres 
output, so x litres input gives 0.42 litres output. 
Similarly, the amount of water reaching node U 
from node B is 0.3y litres. Adding together the 
contributions from the two input nodes gives 
0.42 + 0.3y litres of water output at node U. 


Similarly, we add together the contributions 
from the two input nodes to obtain: 


0.22 + 0.1y litres of water output at node V; 


0.42 + 0.6y litres of water output at node W. 
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Solution 1.3 


(a) 


Inputting 1 litre of water at node U gives 
0.2 litres of water output at node C and 
0.8 litres of water output at node D. 


Inputting 2 litres of water at node V gives 
2 x 0.6 litres of water output at node C and 
2 x 0.4 litres of water output at node D. 


So, in total, we have 0.2 + 2 x 0.6 = 1.4 litres of 
water output at node C,, and 0.8+ 2 x 0.4 
= 1.6 litres of water output at node D. 


Using a method similar to that of part (a) gives 
0.2% + 0.6y + 0.32 litres of water output at 
node C' and 0.82 + 0.4y + 0.7z litres of water 
output at node D. 


Solution 1.4 


(a) 


There are three routes through the network that 
start at node A and end at node D, namely the 
routes through nodes U, V and W. The method 
of calculating the output at node D is to add 
together the contributions from each of the 
possible routes. 


First consider the route that goes via node U. 
From Figure 1.3, 0.4 litres of water flow from 
node A to node U. Also from the figure, if 1 litre 
were input at node U, then 0.8 litres would be 
output at node D. So if 0.4 litres are input at 
node U, then 0.4 x 0.8 = 0.32 litres are output 
at node D. 


Similarly, the amount of water flowing from 
node A to node D via node V is 

0.2 x 0.4 = 0.08 litres, and the amount of water 
flowing from node A to node D via node W is 
0.4 x 0.7 = 0.28 litres. 


So, in total, the amount of water output at 
node D is 


0.32 + 0.08 + 0.28 = 0.68 litres. 


As in part (a), the amounts of water output can 
be found by calculating the flow through each of 
the possible routes. 1 litre of water input at 
node B gives: 

(0.3 x 0.2) + (0.1 x 0.6) + (0.6 x 0.3) = 0.3 litres 
output at node C; 
(0.3 x 0.8) + (0.1 x 0.4) + (0.6 x 0.7) = 0.7 litres 
output at node D. 


Solution 1.5 


(a) The answer can be read off from Figure 1.4: 
0.3 litres of water are output at node C, and 
0.7 litres of water are output at node D. 
(This is much easier than calculating the 
outputs using the combined network, as was 
done in Activity 1.4(b).) 


(b) From Figure 1.4, we find the outputs from each 
input and then add together the contributions. 
This gives 0.32x + 0.3y litres of water output at 
node C, and 0.682 + 0.7y litres of water output 
at node D. 


Solution 1.6 


(a) The first component of the vector is the number 
of litres of water input at node A, which in this 
case is 2. The second component of the vector is 
the number of litres of water input at node B, 
which in this case is 3. So the required vector is 


2 
3}: 
(b) The first component of the vector is 2. 


(c) The three-component vector given in the 
question represents outputs of 0.3 litres of water 
at node U, 0.1 litres of water at node V, and 
0.6 litres of water at node W. 


Solution 1.7 


(a) The vector which represents an input of 1 litre 
of water at node V and no other inputs is 


0 
1 
0 


(Note that the input vector has three 
components because there are three input nodes 
in the network diagram.) 


(b) The input vector has three components. So, 
according to the convention, the matrix 
representing the network must have three 
columns (one for each input node). 


(c) The output vector has two components, since 
there are two output nodes. So, according to the 
convention, the matrix representing the network 
must have two rows. 


(d) By parts (b) and (c), the matrix representing 
the network in Figure 1.6 must have two rows 
and three columns. The number that labels the 
pipe from the first input node (U) to the first 
output node (C) is 0.2. This is the number that 
should go in the first row and first column. The 
number that labels the pipe from the second 
input node (V) to the first output node (C) 
is 0.6, and this should be put in the first row 
and second column. 


SOLUTIONS TO ACTIVITIES 


Proceeding similarly for all the inputs to and 
outputs from the network gives the following 
matrix to represent the network: 


0.2 0.6 0.3 
0.8 0.4 0.7) 7 


Solution 1.8 


1 
(a) The required input vector is | 0 
0 


(b) Start by writing down the matrix representing 
the network and the input vector: 


& 0.6 - 
0.8 0.4 0.7 

0 
Now use the rule which was described in the 
text preceding the activity: for each row of the 
matrix, add together the results of multiplying 
each element in the row with the corresponding 
element in the input vector. In this case, this 
procedure gives 


0.2x1+06x0+03x0\  /0.2 
0.8x1+04x0+0.7x0/  \08/)° 
This is the output vector, which represents an 


output of 0.2 litres of water at node C' and 
0.8 litres of water at node D. 


Solution 1.9 
(a) 0.33 0.25 1\ — (0.383 + 0.25 
0.67 0.75 1) \ 0.67+ 0.75 
_ (0.58 
~ \ 1.42 
(b) 0.33 0.25 - _ (0.33 x 10+ 0.25 x 14 
0.67 0.75 ~ \ 0.67 x 10+ 0.75 x 14 


es 


6.7 + 10.5 


-( 

= (172) 
© (Ger ome) (‘e) -(osrnietonexe) 

( 

( 


12.06 + 4.5 


7.44 
16.56 


i) 
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CHAPTER B2 MODELLING WITH MATRICES 


Solution 1.10 


(a) (i) The first task is to work out which nodes are 
involved in the top right element of the product 
matrix. The top row corresponds to the first 
output, so we are looking for routes which end 
at node C. The right-hand column corresponds 
to the second input, so we are looking for routes 
which start at node B. The network diagram 
with routes from node B to node C highlighted 
is given in Figure $.1. 


Figure S.1 


(ii) The relevant elements in the two matrices 
are the top row of the first matrix and the 
right-hand column of the second matrix, as 
highlighted in Figure S.2. 


O4 0.3 
br 0.6 “a on 


0.1 
08 G4 0.7 04 06 


Figure S.2 


(iii) Adding the products of corresponding 
elements gives 


0.2 x 0.3+ 0.6 x 0.14 0.3 x 0.6 = 0.3. 


So the top right element of the product matrix 
is 0.3, which agrees with the value shown in 
Figure 1.9. 


(b) (i) This time, the bottom right element of the 
product is required; this corresponds to routes 
between the second input node (B) and the 
second output node (D). The network diagram 
with routes from node B to node D highlighted 
is given in Figure $.3. 
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Figure S.3 


(ii) The relevant elements in the two matrices 
are the bottom row of the first matrix and the 
right-hand column of the second matrix, as 
highlighted in Figure S.4. 


08 04 0.7 2 


(i 0.6 a) o¢ & 
04 0.6 


Figure S.4 


(iii) Adding the products of corresponding 
elements gives 


0.8 x 0.3+0.4 x 0.1+0.7 x 0.6 = 0.7. 


So the bottom right element of the product 
matrix is 0.7, which agrees with the value shown 
in Figure 1.9. 


Solution 2.1 


In order to use the method described in the text, the 
sizes of the matrices are written underneath the 
matrices, and the middle two numbers are boxed. 


2 9 oS. 
| hae 
(a) a eer) 
6 8 


3x[22]x3 


The two numbers in the box are equal, so the 
matrices can be multiplied. Ignoring the 
numbers in the box gives the size of the product 
matrix as 3 x 3. (So the product matrix will 
have 3 rows and 3 columns.) 


9 
(b) (1 2 3)]{8 
7 


1x[3 3]x1 


The two numbers in the box are equal, so the 
matrices can be multiplied. Ignoring the 
numbers in the box gives the size of the product 
matrix as 1 x 1. (So the product matrix will 
have 1 row and 1 column. By convention, we 
usually omit the parentheses around a 1 x 1 
matrix and treat it as a number.) 


1 9\ / 0 10 
(c) 7 4 4 4 
—-6 8 —2 7 

sx fae 


The two numbers in the box do not match, so 
these two matrices cannot be multiplied 


together. 
9 

(d) 8 (1 2 3) 
v4 

3x1 1)x3 


The two numbers in the box are equal, so the 
matrices can be multiplied. Ignoring the 
numbers in the box gives the size of the product 
matrix as 3 x 3. (So the product matrix will 
have 3 rows and 3 columns.) 


1 2 3 ne : 
(e) 761 1 0 -1 
2 1 3 

2x[Ea]x:3 


The two numbers in the box are equal, so the 
matrices can be multiplied. Ignoring the 
numbers in the box gives the size of the product 
matrix as 2 x 3. (So the product matrix will 
have 2 rows and 3 columns.) 


Solution 2.2 
2. 4 3.5 11 
(a) AB= {4 2 é eS :) ={/6 10 22 
6 1 t. AF “3h 
3x/2 2/x3 gives 3x3 


2x[3 3]x2 gives 2x2 
2 1\/4 6 10 11 
a € a G $)= (Cs i) 
2x[2 2]x2 gives 2x2 
4 €6\/(2 1 26 4 
pa=(2 9) (3 o)=(2 2) 


2x[2_2|x2 gives 2x2 


“[2. 2) x2 gives 3x2 


SOLUTIONS TO ACTIVITIES 


BA cannot be formed, since the number of 
columns (2) in B does not match the number of 
rows (3) in A. Writing out the sizes of the 
matrices gives 


2x|/2 3/x 2, 


which also confirms that the product BA cannot 
be formed, because the numbers in the box do 
not match. 


Solution 2.3 


(a) The matrix product AB can be formed only if 
the number of columns of A is equal to the 
number of rows of B. Since the number of 
columns (3) of M is not equal to the number of 
rows (2) of M, the product MM cannot be 
formed. 


(b) A must have the same number of rows as 
columns in order for it to be possible to 
calculate its square. 


Solution 2.4 


(a) a= (5 .) and B= (5 _})s then 


5 -l 
2 0 4 101 
(b) a= ( 1 6 5) and B= ( 1 ae 
then 
3.0 5 
A+B=B+A=(_j - :) 
4 0 5 6 2 0 
(c) TA=]6 1 0|]andB=|{|0 -1 4 
2 3 -1 5 0 1 
then 
10 2 5 
A+B=B+A=]6 0 4 
7 3 =0 
Solution 2.5 
are 0.83) (0.36) (1.19 
0.17 0.64)  \ 0.81 
—5 —3 —8 
(ii) 7)+ 1)]= 8 
—A 4 0 


a (1)+(a)+()=() 


(b) (i) Any pair of vectors with different numbers of 
1 
components will do. For example, | 2 | and 
3 


5 cannot be added together, since the first 


vector has three components and the second has 
only two. 
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CHAPTER B2 MODELLING WITH MATRICES 


Solution 2.6 


w (3 B)=-(5 8) 


In each of the answers in this part, the scalar 
multiple is not unique. For example, in part (i), 


ao 8) a (Is 1 
2 15) 4\8 6/)° 
. 1 0 2 0 
(c) since B= (5 1) 2B = (j ) )- Hence 


aep)=(3 4) (4 2)=(2 a): 


Thus A(2B) = 2(AB), as required. 


Solution 2.7 
We obtain 


s-(2 )(3)-G) 


162 C=C), 
fea (41) + a = fae 
“asen()+(1)-() 


A(B+C) = é 7 (3) 7 (G2), 


Thus AB + AC = A(B + C), as required. 
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Solution 2.8 


For the matrix 


(a) 
(b) 
(c) 


P25 i 
B=] 1 i 


bi = it 
ba1 = —4, 
b13 = 0; 


Solution 3.1 


(a) 


S 


The birth rates for juveniles and elderly are zero 
(to 4 decimal places), so we have to consider 
only the number of workers’ babies. 

Rearranging equation (3.1) gives, in millions, 
number of births = birth rate x number of workers 
= 0.0213 x 37.50 

~ 0.7988. 


These data lead to an estimate of approximately 
800 000 births between June 1990 and June 1991. 


To compare the estimated numbers of deaths, we 
multiply the death rate in each category by the 
number in the category. This gives, in millions, 


number of worker deaths = 0.0031 x 37.50 
~ 0.1163, 
number of elderly deaths = 0.0575 x 8.99 
~ 0.5169. 
So there were nearly five times as many elderly 


deaths as those for workers (even though there 
were about four times as many workers). 


Solution 3.2 


(a) 


Table 3.3 tells us that the number that labels 
the pipe from Ap to J; is the proportionate 
birth rate for adults. Table 3.2 gives this a value 
of 0.0172. 


Table 3.3 tells us that the number that labels 
the pipe from Ag to A; corresponds to the 
proportion of adults who do not die in year 0. 
Table 3.2 gives the death rate for adults as 
0.0136, so the proportion surviving is 


1 — 0.0136 = 0.9864 (to 4 d.p.). 


Table 3.3 tells us that the number that labels 
the pipe from Jo to A; corresponds to the 
proportion of surviving juveniles who reach 15. 
Table 3.2 gives the death rate for juveniles as 
0.0008, so the proportion surviving is 

1 — 0.0008 = 0.9992. Using the assumption that 
+ of these juveniles become adults each year 


15 
gives the label as 


75 X 0.9992 = 0.0666 (to 4 d.p.). 


Solution 3.3 


Birth and death rates constant This 
assumption was needed to create the matrix model in 
equation (3.2). It is likely to be valid for only short 
time periods. It is this assumption that limits the 
period over which the model can predict with any 
accuracy. 


No juveniles give birth This assumption was not 
necessary to obtain the type of model given in the 
text. The data in Table 3.2 show that this is a fair 
assumption to make, since the birth rate for juveniles 
is insignificant (0.0000 to 4 d.p.). 


— of surviving juveniles become adults each 
year This assumption was needed to create the 
matrix model. It means that the spread of ages 
within the juvenile subpopulation is uniform. This 
will be the case if the changes in the population are 
small compared with the size of the population. 

A quick check on the predictions of the model will 
reveal the extent to which this is valid. 


No emigration or immigration This factor was 
completely ignored in the text. More data are needed 
to decide whether this is a valid assumption. (In fact, 
this is a fair assumption for the UK population, as 
the rate of immigration is almost equal to the rate of 
emigration.) 


Solution 3.4 
The data in Table 3.2 give the initial values 


Jo\ _ ( 10.92 
Ag) \ 46.49] ° 
Now use equation (3.2) to find estimates for year 1: 
J, \ — (0.9326 0.0172 10.92 \ _ / 10.98 
A, / ~~ \ 0.0666 0.9864 46.49 } ~ \ 46.59 }° 
These estimates for year 1 can be used to calculate 
estimates for year 2: 


Jz\ — (0.9326 0.0172 10.98 \ _ / 11.04 
Az} \ 0.0666 0.9864 46.59 ]} ~~ \ 46.68 }° 
These results can be summarised and interpreted as 


follows. Each of the predictions is expressed in 
millions of people. 


J, ~ 10.98: prediction for juveniles in June 1991. 
A, ~ 46.59: prediction for adults in June 1991. 
Jy ~ 11.04: prediction for juveniles in June 1992. 
A» ~ 46.68: prediction for adults in June 1992. 


(The values of Jz and Ag were calculated using the 
values of J; and A; produced by the calculator before 
these were rounded. If you used rounded values, then 
your answers may be slightly different. The more 
years calculated in this way, the bigger the 
difference.) 


SOLUTIONS TO ACTIVITIES 


Solution 5.1 
The products AI, IA and IB can be formed, with 


the following results: 
3°44 
(4 t)-4 


u=($ 16%) 
GL a 4 
eC 401 ee 


The product BI cannot be formed (because B has 
3 columns and I has 2 rows). 


Solution 5.2 


(a) According to the assertion, Ib = b since the 
product Ib can be formed (I has 3 columns and 
b has 3 rows). 


(b) Verification: 


1 0 0 by 
b=} 0.1 0} |b 
OO AF \bs 
1 x by 0 x be 0 x bg 


= 0 x b) +1 x bg +0 x b3 
0x 6b) +0 x bg +1 x bg 


Solution 5.3 


‘a 
Q 
iw) 
| 

a 
ee © 
| 
ar 
a ™~ 
| 
Re bh 
or 
Ne oe"  ~ 
lI 
Pa ¥ 
oF 
He © 
So 


(iii) Pa=(j 7 é : =(4 :)) 
ap=(4 ») (4 a) =(c ») 


(b) Each pair of matrices in parts (a)(i) and (a) (ii), 
taken in either order, gives the identity matrix 
as product. 


The pair of matrices in part (a)(iii), taken in 
either order, gives a product which looks rather 
like the identity matrix. In fact, the product 
matrix is a scalar multiple of the identity 
matrix, namely 3I. 
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CHAPTER B2 MODELLING WITH MATRICES 


Solution 5.4 
Multiplying the two matrices together gives 
a bd d —b\ _ fad—bc 0 
c d -c aj 0 —bc+ ad }° 


So the product is (ad — bc)I, the 2 x 2 identity 
matrix scalar multiplied by ad — bc. 


Solution 5.5 
(a) det A = 26 — 25 =1#0; hence A is invertible 


and 
“1 2 -5 
7 = 13 )° 


(b) detB =6—0=6 ¥ 0; hence B is invertible and 


2 0 3 (0 
i _ {3 
oei(? = (1°) 
6 2 
(c) det C = 1-—1=0; hence C is non-invertible, 


that is, C~! does not exist. 


(d) det D = 2.25 — 1.25 = 1# 0; hence D is 
invertible and 


Solution 5.6 


(a) To write the equations 


2x + 3y = 3, 
xr+4y=-1, 
in matrix form, first write down the coefficient 
matrix ; j . Then write down the matrix 


form of the equations: 


(ra) G)=Ca) 


(b) The matrix form of the equations 


22 =6, 
3x + 6y = 15, 


(6) (3) = (is) 


(c) The matrix form of the equations 


is 


3 
5 
4 a 

gor gy——L 


a— 2y = 18, 


is 


Pr 
ow Oe 


ye 


Solution 5.7 


(a) Activity 5.6(a) gives the matrix form of the 
equations as 


2. 3 oo 3 
1 4 y) \-1)° 
The determinant of the coefficient matrix is 


2x4-3x1=5, 


so the inverse of the coefficient matrix is 


The solution is z = 3, y = —1. (You should 
always check your solution by substituting back 
into the original equations. ) 


(b) Activity 5.6(b) gives the matrix form of the 
equations as 


2 0 x\ (6 
36 y) \15)° 
The determinant of the coefficient matrix is 
2x6-0x3=12, 


so the inverse of the coefficient matrix is 


cf @ O\. 7 go 8 
m\-3 a} 2 a] 
Hence 
6 
15 


6+0x 15 
x6+ 5x15 


ale © 


NP RIF ble 


II 
fn 
“ar 
AIR 
“—" 4 


The solution is = 3, y= 1. 


(c) Activity 5.6(c) gives the matrix form of the 
equations as 


G)G)-). 


ou ooo 


The determinant of the coefficient matrix is 


5x 5—(-3)xs=1 


so the inverse of the coefficient matrix is 


oo oe 


lI 
mo Ne OR eS 
| —— 


reo Negi oe OU oe 
ane ANS oy 7 


The solution i 


mn 


z=10,y=-15. 


Solution 5.8 


(a) The matrix form of the equations is 


2 -3 x\ (5 
-4 6 y) \T7)° 
The determinant of the coefficient matrix is 
2 x 6 — (—3) x (—4) = 0, 


so the matrix has no inverse. Hence the matrix 
method does not provide a solution. 


Now let us try solving the equations using the 
elimination method. Start by doubling the first 
equation in preparation to eliminate z: 
4x — 6y = 10, 

—4a + 6y = 7. 
Now add the two equations to eliminate x; this 
gives 

Ox — Oy = 17. 
Obviously, the equation 0 = 17 has no solution, 


so the original equations also have no solution. 


The matrix form of the equations is 
1 2 o\ (£6 
—3 -6 y) \ 18)° 
The determinant of the coefficient matrix is 


1 x (—6) — 2x (-3) =0, 


so the matrix has no inverse. Hence the matrix 
method does not provide a solution. 


SOLUTIONS TO ACTIVITIES 


The elimination method can start by tripling the 
first equation to give 


3x2 + 6y = —-18, 
—3x — by = 18. 


Now add the two equations to eliminate x; this 
gives Ox + Oy = 0, that is, 0 = 0. So the 
elimination method fails to find values for x 
and y, and this is the same outcome as with the 
matrix method. 


(In fact, the equation 0 = 0 is true for all values 
of x and y, and this is an indication that the 
equations have infinitely many solutions, as 
pictured in Figure 5.3, following the activity.) 


Solution 5.9 


(a) 


The matrix form of the equations is 


(1 4) G)-@) 

-1 4 y 7). 

The determinant of the coefficient matrix is 
5 x 4— (3) x (-1) = 17. 

This is non-zero, so the coefficient matrix is 


invertible and the equations have a unique 
solution. 


The matrix form of the equations is 
(5 3) (4)= Gs) 
-9 3 y 18 } © 
The determinant of the coefficient matrix is 
3x 3—(-1) x (-9) =0. 


This is zero, so the coefficient matrix is 
non-invertible and the equations do not have a 
unique solution. (The equations represent the 
same line.) 


The matrix form of the equations is 


6 —2 u\ 1 
3-7 yJ) \-1}¢ 
The determinant of the coefficient matrix is 
6 x (—7) — (—2) x 3 = —36. 


This is non-zero, so the coefficient matrix is 
invertible and the equations have a unique 
solution. 


De) 


Solutions to Exercises 


Solution 1.1 
. (0.35 0.85) (2 
ay) & i. (5) 
_ (0.35 x 2+ 0.85 x 6 
0.65 x 2+0.15 x 6 


_ (07454 
ae err ain: 
_ (58 
=e 
iy (935 0.85) ( 50 
0.65 0.15) \ 100 


0.35 x 50 + 0.85 x 100 
0.65 x 50 + 0.15 x 100 


17.5 + 85 
32.5+ 15 
102.5 
47.5 
(iii) 0.35 0.85 60 
0.65 0.15 40 


_ {0.35 x 60 + 0.85 x 40 
~ \ 0.65 x 60+ 0.15 x 40 


_ (21434 
~ \ 3946 


=() 


(b) The network is shown in Figure S.5. 


Figure S.5 


Solution 2.1 


(a) All the products can be formed except CA 
and D?. 


8 12 10 17 8 
(b) AB= (1 9 a cp=({ 


4 ae Sai 
DA=]|8 , DB=|{ 16 24 20], 
3 6. 2-2 


54 


Solution 2.2 


(a) All the combinations can be formed except 
A+BandC+A. 


0.5 0.3 0.2 
(b) | 01}]-—[ 04] =| -03 
0.4 0.3 0.1 


Solution 3.1 


Let A, and J,, denote the numbers (in millions) of 
adults and juveniles, respectively, at the start of 
year n, where n = 0 corresponds to the year starting 
in January 2000. 


(a) The same reasoning as that used for the UK 
population leads to the following calculations of 
numbers labelling pipes. 


Pipe 
From To Label 
Jo Ji F4(1 — 0.0139) ~ 0.9204 
Jo Ay (1 -0.0139) ~ 0.0657 
Ap Ji 0.0520 
Ap Ay 1 — 0.0231 ~ 0.9769 


These labels have been added to the network 
diagram shown in Figure $.6. 


do Ao 


0.9204 ¥ Y 0.9769 


dy Ay 


Figure S.6 


(b) In41 \ _ (0.9204 0.0520 In 
An+i)  \0.0657 0.9769 An, 


(n =0,1,2,...) 


Solution 5.1 
(a) det A = (—3) x 4—(-5) x 1 = -12+5 = -7, 


so 


so 


gaut/2—4)\.f 9s 
Ge ee ee 


so 


3 4 3 

cnt ae 

2 4 3 4 
1\-$ § = 


(Note that C-! = C.) 


cnlos US 
Se 


Solution 5.2 


(a) The matrix form of the equations is 


(1 4) G)-@)- 


From Exercise 5.1(a), the inverse of ¢ 4) 


_4 _5 

F 7 7 
is 1 3}? 80 

7 7 


The solution is = —2, y= 1. 


SOLUTIONS TO EXERCISES 


(b) The matrix form of the equations is 


(3 -1)@)= (io) 


The inverse of A = é 7 is 


ate f-1 -\_fm 3 
-12\-3 6 , 2h 
sO 
c\ 3 é 4 
vy} ~\ a 4] \10 
( on 
qx 4+ (—4) x10 
iL 5 
| ers 
1-5 
_{ 2 
=(_4]- 
The solution is « = 2, y = —4. 


(c) The matrix form of the equations is 


(2)()-(2) 


The inverse of A = (1 a is 
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coefficient matrix 40 
column of a matrix 11 
components of a vector 10 


determinant 38 
determinant test 38 
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element of a matrix 11 
equality of matrices 18 


identity matrix 35 
inverse of a2 x 2 matrix 36, 37 
invertible matrix 38 


leading diagonal of a matrix 36 


main diagonal of a matrix 36 
matrix 11 

mxn 18 
matrix addition 21, 25 
matrix multiplication 16, 19, 25 
matrix subtraction 23 
matrix—vector multiplication 13 


network, matrix representation of 10 
network diagram 6 

node of a network 6 

non-invertible matrix 38 


pipe of a network 6 
power of a matrix 21 
product matrix 16 


row of amatrix 11 


scalar 23 

scalar multiplication 23 
size of a matrix 18 
square matrix 21 


vector 10 
vector addition 22 
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